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Abstract 

The charge exchange of vector polarised deuterons on a polarised hydrogen target has been studied in a high statistics experiment 
at the COSY-ANKE facility at a deuteron beam energy of Td = 726 MeV. By selecting two fast protons at low relative energy E pp , 
the measured analysing powers and spin correlations are sensitive to interference terms between specific neutron-proton charge- 
exchange amplitudes at a neutron kinetic energy of T n ss |Td = 363 MeV. An impulse approximation calculation, which takes into 
account corrections due to the angular distribution in the diproton, describes reasonably the dependence of the data on both E pp 
and the momentum transfer. This lends broad support to the current neutron-proton partial-wave solution that was used in the 
estimation. 
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It is a consequence of the nucleon spins that, assum¬ 
ing charge independence, five complex amplitudes are 
needed to describe neutron-proton elastic scattering [1]. 
This means that, above the pion production threshold, 
at least nine independent measurements are required at 
each scattering angle to allow an unambiguous partial 
wave decomposition. Some of the resulting observables, 
which could depend on up to three spin projections [2], 
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are difficult to determine and values may only be obtained 
indirectly through combinations of other measurements. 

It was shown several years ago [3] that, at small momen¬ 
tum transfers between the deuteron and the diproton, the 
tensor analysing power in the deuteron charge exchange on 
hydrogen, dp -A { pp} s n , is closely linked to the spin trans¬ 
fer in neutron-proton large angle scattering, pn —> ftp, pro¬ 
vided that the excitation energy E pp in the final diproton 
is very low. Due to the Pauli principle the two protons are 
then dominantly in the 1 Sq state with antiparallel spins, 
here denoted by {pp} s , so that there is then a spin-isospin 





flip to this state from the initial deuteron. 

Further information on the neutron-proton scattering 
amplitudes can be obtained through measurements of the 
analysing powers and spin correlations in the dp —> {pp} s n 
reaction and measurements of this type were carried out 
at deuteron beam energies of Td = 1.2 GeV and 2.27 GeV 
to investigate the np amplitudes at neutron kinetic ener¬ 
gies of T n ft ^Td = 600 and 1135 MeV [4]. Both transverse 
spin correlations and the proton and deuteron analysing 
powers were investigated and the results were found to be 
consistent with modern partial wave solutions [5] at T n = 
600 MeV, while failing badly at 1135 MeV. We report here 
on a similar investigation carried out at Td = 726 MeV in 
a high statistics experiment, where tighter cuts could be 
placed on E pp and small effects could be studied in detail. 
Since one might expect that the np partial wave amplitudes 
should be fairly reliable at 363 MeV, this is the ideal test¬ 
ing ground to establish quantitatively the validity of the 
theoretical modelling of deuteron charge exchange [6]. 

The experiment was undertaken using the ANKE mag¬ 
netic spectrometer installed at an internal target position 
of the Cooler Synchrotron (COSY) at the Forschungszen- 
trum Jiilich [7]. Data were taken in parallel with those 
used to determine the spin correlations in quasi-free np —> 
{pp}sir~ [8] and more details of the experimental proce¬ 
dure, in particular of the measurements of the beam and 
target polarisations, are to be found in this reference. 

Only deuteron beams with vector polarisation were used 
in this experiment and these had ideal values of p\ = | 
and p^ = — |. The polarisations measured at the injection 
energy of 75.6 MeV with the low energy polarimeter were 
p\ = +61 ± 4% and p ^ = — 50 ± 3% for the two states, while 
the tensor polarisations were shown to be below 2%. 

In order to increase the luminosity in the experiment, a 
jet of polarised atomic hydrogen was fed into a 25 p thick 
teflon-coated aluminum storage cell target with dimensions 
x x y x z = 15 X 19 x 390 mm 3 . Here the y-direction is 
perpendicular to the COSY plane and the ^-direction is in 
this plane but perpendicular to the beam (z) direction. The 
polarisation, p p , of the target was in the y direction and 
its sign was reversed every five seconds. The mean value 
of the polarisation was determined through the study of 
the quasi-free np —> dir 0 asymmetry to be \p p \ = 69% ± 
2%(stat) ± 3.5% (syst). A more precise value of the product 
of the magnitudes of the beam and target polarisation was, 
however, extracted from an analysis of the pion production 
data themselves, which gave an averageoflpdUppI = 0.373± 
0.015 [8] % 

Although the ANKE spectrometer is equipped with other 
elements, the only detector used in the charge-exchange ex¬ 
periment was the forward detector (FD) that identified and 
measured the two fast final protons from the dp {pp}n 
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reaction or, for polarisation studies, the fast deuteron and 
spectator proton from the dp —> dpir 0 reaction. The FD 
comprises a set of multiwire proportional and drift cham¬ 
bers and a two-plane scintillation liodoscope [9]. 

Having registered two charged particles in the FD, the 
isolation of the dp —> ppn reaction depends on identifying 
these as protons on the basis of time-of-flight criteria that 
are described in detail in Ref. [4]. For this purpose the 
difference in the times of flight of the particles recorded in 
the FD is compared to that calculated on the assumption 
that the two particles are both protons. This procedure 
suppresses enormously the background that is associated, 
for example, with deuteron-proton pairs coming from dp 
elastic scattering. 
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Fig. 1. Distribution in missing masses for the dp —»• ppX reaction for 
726 MeV deuterons incident on a storage cell filled with polarised 
hydrogen. Events that originated from the walls of the cell gave rise 
to the tail at high Mx • The shape of this background was simulated 
by filling the cell with nitrogen gas and this led to the shaded area. 
The solid line represents the fit of a Gaussian plus the numerical 
values of the scaled background histogram. 

The missing-mass Mx distribution of the identified dp — > 
ppX reaction shows a striking peak around the mass of the 
missing neutron, as illustrated in Fig. 1. The long tail to 
higher missing masses arises from events originating from 
the walls of the storage cell. The shape of this background 
was simulated by filling the cell with nitrogen gas. This gave 
rise to the shaded area in the figure which, after normalising 
the distribution at high Mx , could be reliably subtracted 
bin by bin. The resulting dp —> ppn events were then placed 
in 20 MeV/c bins in the momentum transfer q between 
the deuteron and diproton and 2 MeV bins in the diproton 
excitation energy E pp . 

For a vector polarised deuteron beam incident on a po¬ 
larised hydrogen target, where both polarisations are in the 
y direction, the ratio of the numbers of polarised N(q, <f>) 
to unpolarised N°(q) events has the form [10]: 

= 1 + Pp a I (<?) cos 0 + | Pd,Ay (q) cos (j> + 

+ |my[(1 + cos2(j)) Cy,y{q) + (1 - cos2 <j>) C XtX (q)\, (1) 
where the azimuthal angle ( j> is measured from the axaxis. 
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Fig. 2. Values of the deuteron (vector) analysing power Ay (black in¬ 
verted triangles) and proton analysing power Ay (blue open circles) 
for the dp — t {pp}n reaction at 726 MeV. Data were placed in bins of 
(a) E pp < 2 MeV, (b) 4 < E pp < 6 MeV, and (c) 8 < E pp < 10 MeV; 
only statistical errors are shown. The curves are impulse approxima¬ 
tion estimates [6] folded with experimental acceptance. These used 
the current SAID neutron-proton partial wave solution [5] as input. 
The dashed lines neglect the predicted dependence on the angle be¬ 
tween the relative momentum in the diproton and the momentum 
transfer. The model predictions for this are included in the solid lines. 

From studying the <f> dependence of the count rates for 
the four combinations of beam and target polarisations it is 
possible to extract separately the values of the proton and 
deuteron vector analysing powers as well as the two spin 
correlations. The results for these observables are shown in 
Figs. 2 and 3. The deuteron vector analysing power Ay of 
Fig. 2 remains very small over our whole q range, only (pos¬ 
sibly) exceeding 1% in magnitude for q > 120 MeV/c. The 
proton analysing power A?, though small, is much larger 
than Ay. These two features are very similar to the results 
found at 600 MeV per nucleon [4], though the statistical 
precision of the current data is much higher. 

Due to the ANKE exit window being much wider in the 
horizontal direction than in the vertical, for the larger val¬ 
ues of q the data are more populated near sin0 = 0. It 
follows from Eq.(l) that, in this limit, the spin correlation 
Cy t y is better measured than C XiX , and this is seen in Fig. 3. 
However, in order to assess the significance of these results 
we must turn to a reaction model. 

In impulse approximation the amplitude for the dp —> 
{pp}n charge-exchange reaction is proportional to the 
np —> pn charge-exchange amplitude times a form factor 


Fig. 3. Values of the spin correlations O y , v (black inverted trian¬ 
gles) and Cx.-.r (blue open circles) for the dp {pp}. s n reaction at 
726 MeV. The conventions are identical to those of Fig. 2. 

that represents the overlap of the initial deuteron wave 
function with that of the outgoing diproton [3,6]. The el¬ 
ementary np —> pn amplitude may be written in terms of 
five scalar amplitudes in the np c.m. frame as: 

fnp = a(q) + ij(q)(cri + er 2 ) • n + P{q)(cri ■ n )(er 2 • n) 

+<5(q)(cr 1 • m)(<r 2 • m ) + e(q)(cr 1 ■ l)(cr 2 ■ l), (2) 

where q = y/—t is the three-momentum transfer between 
the initial neutron and final proton and the Pauli matrices 
u are sandwiched between neutron and proton spinors. Of 
the unit basis vectors, l lies along the mean of the initial 
proton and final neutron momenta, n lies along q , and m = 
nxl.lt should be noted that the amplitudes of Eq. (2) are 
actually linear combinations of the standard elastic ones, 
defined for example in Ref. [2]. This important distinction 
arises because the spin dependence that is made explicit 
here is that corresponding to charge exchange [3,11]. 

Values of the amplitudes of Eq. (2) at 363 MeV can be ex¬ 
tracted from the current partial wave solution of the SAID 
group [5] and the ones that are relevant to the current work 
are shown in Fig. 4 as functions of the momentum transfer q. 
Since only relative phases are significant in the discussion, 
these amplitudes have been rotated in the complex plane 
to make /? real for all q. Apart from the obvious features 
that f3 = 5 and 7 = 0 at q = 0, the most notable behaviour 
is the zero in the <5 amplitude at q k, (140 — lOi) MeV/c. 
This amplitude is strongly influenced by one pion exchange 


Q 





























0 100 200 


q [MeV/c] 


Fig. 4. Values of the (a) real and (b) imaginary parts of the np 
amplitudes at 363 MeV predicted on the basis of the current SAID 
PWA solution [5]. The amplitudes of Eq. (2), normalised to da / dq 2 , 
have been rotated such that /3(q) is real. The curves are for /3 (black), 
5 (blue), e (red), and 7 (magenta). It should be noted that the scales 
in panels (a) and (b) are very different. 


and, in the simplest distorted model, this has a zero when 
q = m^c, where m w is the mass of the charged pion. 

Although the resulting dp —> {pp}n amplitudes will be 
evaluated taking into account higher partial waves in the 
pp system, using an update of the program of Ref. [ 6 ], it 
is useful for a qualitative discussion to consider the results 
that follow if one retains only the pp hS’o configuration that 
dominates at low E pp . In this case the impulse approxima¬ 
tion model predicts that: 


0 , 

A? = —2/m(/3*7)/(l/?| 2 + ItI 2 + M 2 + l<5| 2 ), 

C v , y = -2Re(e*6)/m 2 + h\ 2 + \e\ 2 + \5\ 2 ), 

C x , x = -2Re(e*P)/m 2 + 111 2 + k| 2 + l<5| 2 ). (3) 


In addition to taking the higher pp partial waves into ac¬ 
count, away from q = 0 these formulae have to be mod¬ 
ified to include the effects of the deuteron D -wave and 
the Wigner rotation that arises from a change in reference 
frame [3]. 

In the 1 <S'o limit of the impulse approximation the spin 
correlations are linked at q = 0 to a combination of neutron- 
proton spin-correlation and spin-transfer parameters, as 
defined in Ref. [2], through 


C'x.x(O) = C„,„(0) 


2[AoOnn( 7T ’) -^roOno(n)] 

3 Aozzo(n) 2ATQ nn o(if) 


( 4 ) 


The curves shown in Figs. 2a and 3a represent the full 
impulse approximation calculations for E pp < 2 MeV [ 6 ]. 
Though in this case the two protons should be dominantly 
in the kSo configuration, there might still be some small 
deviation from the A p = 0 prediction of Eq. (3). However, 


both the data and the predictions are at the 1 % level and it 
is hard to draw firm conclusions in view of the systematic 
uncertainties. The situation is much clearer for the proton 
analysing power of Fig. 2a, where the current SAID solu¬ 
tion provides a quantitative description of the experimental 
data. Since the real and imaginary parts of 7 in Fig. 4 are 
of comparable size, this suggests that the phase between j3 
and 7 is well reproduced in the partial wave solution [5]. 

The general shapes of the C XtX and C y ^ y predictions in 
Fig. 3a are very much in line with the measured data though 
there are some quantitative differences. C y ^ y changes sign, 
as one would expect from the simple one-pion-exchange 
contribution to the S amplitude, though this happens at a 
few MeV/c higher than the prediction. While being inde¬ 
pendent of the uncertainties in the beam and target polari¬ 
sations, this crossing is a very sensitive test of the S/e inter¬ 
ference because of the presence of the small imaginary part 
in the zero of the S amplitude at q ~ (140 — 10?') MeV/c. 

The prediction of C XjX ( 0) ~ —0.62 at q = 0 is to be 
compared to the value of ~ —0.57 extracted from the ex¬ 
perimental data for E pp < 2 MeV. The significance of the 
discrepancy here must be judged against the systematic 
uncertainties in the experiment and the modeling. System¬ 
atic effects in the data could arise, for example, from the 
choice in fitting limits combined with some imperfection 
in the background description, but these can be estimated 
conservatively to be below 0.01 for the analysing powers 
and 0.03 for the spin correlations. To these must be added 
the 4% associated with the product of the beam and target 
polarisations. The uncertainties arising from the beam or 
target polarisations are small compared to the statistical 
errors in A y and A?. 

The SAID single-energy solution yields error bars on 
the np spin-transfer and correlation parameters needed to 
evaluate Eq. (4). The resulting uncertainty in C XiX ( 0) of 
±0.024 [12] is very much a lower limit because it does not 
include any uncertainties in the model assumed in the SAID 
analysis or in the data selection [5]. It must also be stressed 
that few of the four np observables appearing in Eq. (4) 
have been directly measured near 180°. Furthermore, the 
SAID predictions do show very strong angular dependence, 
which makes any extrapolation in angle less reliable. On 
the other hand, the uncertainty of 4% in the product of the 
beam and target polarisations would correspond to a sys¬ 
tematic error of ±0.025 in the determination of C XiX ( 0). 

The dp —> {pp} s n data were here described using a plane- 
wave impulse approximation [3,6] and the largest correction 
to this picture is likely to come from double scattering inside 
the deuteron. Evaluating the kSo contribution in an eikonal 
approach, it has been shown that the spin dependence of 
the deuteron tensor analysing powers is little changed by 
this correction for q < 140 MeV/c [13]. An estimation of 
the double scattering [3,13] at q = 0 indicates that this 
modifies the prediction for C X)X ( 0) by only 0.003. This is 
an order of magnitude less than the quoted uncertainties 
and so can be safely neglected. 










For larger values of E pp , P and higher waves become 
significant so that the angular distribution in the diproton 
is no longer isotropic. If k is the relative momentum in the 
diproton, then E pp = k 2 /m, where m is the proton mass. 
At large values of k and q, where the effects of the Pauli 
exclusion principle and final state interactions are small, 
there will be a quasi-free peak at k = ±q/2. Quite generally 
therefore, away from the small E pp region there will be a 
significant dependence on the angle 9p q between k and q [3]. 
It is important to note that this effect is already included 
in the computer program of Ref. [6]. 

Panels b and c of Figs. 2 and 3 show the experimental re¬ 
sults obtained in the bins 4 < E pp < 6 MeV and 8 < E pp < 
10 MeV, respectively. The broken lines indicate the plane 
wave impulse approximation prediction where one ignores 
the angular distribution in 9k q ■ When the predictions [6] 
of this angular dependence are included one obtains the 
solid lines in the figures. Though most of the changes are 
small compared to the uncertainties in the np input data, 
these generally go in the right direction, especially for the 
deuteron vector analysing power. The same is also true for 
the results in two E pp bins that are not shown here. 

Spin correlations and analysing powers have been mea¬ 
sured in deuteron charge exchange on hydrogen, dp —> 
{ pp}n , at a beam energy of 726 MeV. The high statistics of 
this experiment allowed tight cuts to be placed on the pp ex¬ 
citation energy E pp . The agreement of these data with the 
impulse approximation model at very low E pp , where the 
% state will be dominant, shows that the np amplitudes 
obtained from partial wave analysis [5] must be broadly 
correct. There are slight discrepancies near the forward di¬ 
rection but these are of such a size that they could originate 
from the partial wave solution or from uncertainties in the 
experimental data presented here. 

At larger values of E pp there are small effects associated 
with P and higher partial waves in the pp system that lead 
to some non-isotropy in the diproton angular distribution. 
These were studied by putting the data in five 2 MeV bins 
in E pp . Though the acceptance in ANKE is less complete at 
large vales of E pp , the plane wave impulse approximation 
model describes all these effects. 


We are grateful to other members of the ANKE Col¬ 
laboration for their help with this experiment and to the 
COSY crew for providing such good working conditions, es¬ 
pecially in respect of the polarised beam. The values of the 
SAID neutron-proton amplitudes were kindly furnished by 
1.1. Strakovsky. This work has been partially supported by 
the Forschungszentrum Jiilich COSY-FFE #73 and #80, 
and the Georgian National Science Foundation. 

References 

[1] M.H. McGregor, M.J. Moravcsik, and H.P. Stapp, Ann. Rev. 
Nucl. Sci. 10 (1960) 291. 

[2] J. Bystricky, F. Lehar, P. Winternitz, J. Phys. (Paris) 39 (1978) 
1 . 

[3] D.V. Bugg, C. Wilkin, Nucl. Phys. A 467 (1987) 575. 

[4] D. Mchedlishvili et al., Eur. Phys. J. A 49 (2013) 49. 

[5] R.A. Arndt, LI. Strakovsky, R.L. Workman, Phys. Rev. C 62 
(2000) 034005; R.A. Arndt, W.J. Briscoe, 1.1. Strakovsky, 

R.L. Workman, Phys. Rev. C 76 (2007) 025209; 

http://gwdac.phys.gwu.edu. 

[6] J. Carbonell, M.B. Barbara, C. Wilkin, Nucl. Phys. A 529 (1991) 
653. 

[7] S. Barsov et ah, Nucl. Instrum. Meth. A 462 (1997) 364. 

[8] S. Dymov et al., Phys. Rev. C 88 (2013) 014001. 

[9] S. Dymov et al., Part. Nucl. Lett. 2 (119) (2004) 40. 

[10] G.G. Ohlsen, Rep. Prog. Phys. 35 (1972) 717. 

[11] F. Lehar, C. Wilkin, Phys. Part. Nuclei Lett. 7 (2010) 235. 

[12] R. Workman, private communication (2014). 

[13] C. Wilkin in NN and 3N Interactions, Ed. L. Blokhintsev and 
I.I. Strakovsky, (Nova Science, N.Y., 2014). 



